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Abstract 



The effect of a random phase screen on laser beam wander in a 
turbulent atmosphere is studied theoretically. The method of photon 
distribution function is used to describe the photon kinetics of both 
weak and strong turbulence. By bringing together analytical and nu- 
merical calculations, we have obtained the variance of beam centroid 
deflections caused by scattering on turbulent eddies. It is shown that 
an artificial distortion of the initial coherence of the radiation can 
be used to decrease the wandering effect. The physical mechanism 
responsible for this reduction and applicability of our approach are 
discussed. 
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1 Introduction 



When a beam of light propagates through the turbulent atmosphere of the 
Earth, it experiences random fluctuations in the refractive index. Fluctua- 
tions of the refractive index are due to turbulent eddies caused by stochastic 
variations of the temperature. The characteristic scales of the atmospheric 
inhomogeneities range from millimeters (the inner radius of the eddies, Iq) 
up to one hundred meters (the outer radius of the eddies, L ). Those inho- 
mogeneities which are large compared with the diameter of the beam tend to 
deflect the beam, whereas those inhomogeneities which are small compared 
with the beam diameter tend to broaden the beam but not deflect it signif- 
icantly. As a result we can observe a broadened laser spot whose centroid 
randomly moves because of the motion of individual eddies. The average 
beam radius is determined by the overall scattering effect, i.e. by both the 
beam broadening and the centroid wandering averaged over a sufficiently 
long time. 

Beam wandering, as well as the scintillation index, is an important char- 
acteristic of the radiation determining its utility for practical applications 
(for example, for purposes of uninterrupted laser tracking and pointing). 
Thus we will study here the possibility of controlling this effect by means of 
artificially decreasing the initial coherence of the radiation using a random 
phase screen. This screen introduces random (spatial and temporal) phase 
distortions into the wave front of the exiting beam. Therefore, after passing 
the phase screen, the initially coherent laser beam becomes partially coher- 
ent. Its coherence length, l c , in the direction perpendicular to the direction 
of propagation becomes smaller than the diameter D of the aperture. As 
a result the initial angular spread of the beam, that is due to diffraction, 
increases from X/D to X/l c , where A is the wavelength of the radiation. (See, 
for example, Refs. p]-[3].) From the viewpoint of possible laser applications, 
the beam broadening is a negative factor that reduces the intensity of radi- 
ation field. At the same time, the wandering effect can become smaller just 
due to the broadening. 

The above comments concern the case of not too long propagation paths 
when the diffraction broadening (which does really depend on the partial co- 
herence) dominates over broadening caused by the atmospheric turbulence. 
But there is another important effect of the phase screen on the statistical 
properties of the radiation propagating in the atmosphere. It is shown in 
papers [I]-[8J that the decrease of the initial coherence may result in lowering 



2 



the normalized variance of the intensity (i.e. the scintillation index) even 
in the case of strong turbulence. This effect takes place only for the case 
of a "slow" detector. The term "slow detector" means the detector has an 
integration time greater than the characteristic time of phase variation intro- 
duced by the phase screen. Since the suppression of the intensity fluctuations 
is of great practical importance, it is also interesting to study the behavior 
of the beam wandering (which is also expressed in terms of local fluctuations 
of the irradiance intensity) for the same experimental arrangement. Thus, in 
what follows the importance of a random phase screen for the case of strong 
turbulence will be elucidated. 

To describe the effect of beam wandering, we will use here an approach 
based on the photon distribution function [B], [5]. 



2 Theoretical description and calculations of 
the wander effect 

The position of beam centroid, H w (z,t), is determined by the expression 

Jdr ± T ± I(T,t) 
* w[Z,) fdr ± (I(r,t))> [i) 

where (...) means averaging over different realizations of the refractive index 
inhomogeneities and source fluctuations, r = {r±,z}, r± = {x,y}, the z-axis 
is along the initial direction of the beam propagation; and the coordinate 
r = corresponds to the center of the exit aperture. 

Following Ref. [6J we express the intensity of photon flux J(r, t) in terms 
of the photon distribution function /(r, q, t) as 

I(r,t) =c$>u> q /(r,q,t), (2) 
q 

where c is the speed of light in a vacuum, u q = cq, 

/(r,q,t) = ^E^ kr ^ +k/2 W, (3) 

k 

V = L x L y L z is the normalizing volume, and 5+ and b n are the creation 
and annihilation operators of photons with the wave-vector q. The operator 
function /(r, q, t) describes the photon density in (r, q) space at time t. For 
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a detailed description of the photon field in a beam with radius Rb, it is 
sufficient to restrict the sum in Eq. ([3]) with some value k , |k| < k Q , where 
R^ 1 « k Q « A -1 . In this case, the distribution function obeys the kinetic 
equation (see Ref. [6]) 

{«9 i + c q d r + F(r)d q }/(r,q,t) = 0, (4) 

where c q = du q /dq, F(r) = u d r n(r), and n(r) is the fluctuating constituent 
of the atmospheric refractive index ((n(r)) = 0, |n(r)| << 1); uo = cqo is the 
central frequency of laser radiation. 

Using Eqs. ([T]) and ([2]), we can easily obtain the variance of R w (z,t) if 
one knows the correlation function of the distribution function (//). The 
analysis is very simple for the case of weak turbulence (or short propagation 
distance). In this case, it is convenient to use not directly Eq. ([1]), but a 
modified expression for it. The following consideration is in the spirit of 
Cook's approach [5] who has used the similarity between a parabolic equa- 
tion describing paraxial optical beam and the Schrodinger equation. The 
application of the Ehrenfest's theorem has made it possible to develop an 
approximate method to study beam wandering effect in [9J. In contrast to 
Cook, we proceed not from Ehrenfest's theorem, but from the definition ([1]) 
and equation (JH). A simple relationship between the beam centroid displace- 
ment and the refractive index fluctuations can be easily obtained within our 
formalism based on the kinetic equation Eq. (jlj) 

(d z + - c d t yR(z,t) = J dr ± J2^f(r, q ,t)(J dv^fi^^t))" 1 (5) 

for the case of a stationary beam. Here and in what follows, the paraxial 
approximation (|q — qo| << qo ) is assumed throughout the beam trajectory. 

In the lowest order with respect to fluctuating refractive index, the de- 
pendence of / on n(r) in ([5]) has to be neglected. Therefore, the variance of 
beam wandering is given by 

(R*) = f f dz ± dz 2 (z - z x )(z - z 2 )x 
Jo Jo 

J d^dr' ± E ( ^ 8 f^ ) (/(r, q, t)f(r>, q', t))( / dr ± £ /(r, q, t))'* , 

(6) 



4 



where r = {rj_, Z\}, r ; = {r'j_, z 2 }, and /(r, q, t) satisfies Eq. (jlj) with F = 0. 
The distribution function at time t can be expressed via its value at the 
instant of photon exit from the source, to, as 

/(r,q,t) = /[r-c q (t-t ),q,t = to] = - £ e-*^^ &J +k/2 6 q _ k/2 | t=to , 

k 

(7) 

where t — to = z/c. In what follows we will put to = for simplicity. 

Two independent averagings should be undertaken in the right-hand part 
of Eq. ([6]). These are averaging over the turbulence configurations, (nn), 
and averaging over the fluctuations of the source, including fluctuations in- 
troduced by the random phase screen, (//). The first of these is determined 
by the known [TD] expression 

(n(r)n(r')) = J dge^-^g), (8) 
where the explicit term for ijj is given by 

</>(g) = 0-033C n [g2 + L -2 ]U/6 ■ (9) 

Eq. ([9]) is referred to as the von Karman spectrum. 

The other averaging accounting for the effect of a "slow" detector is given 
by 0: 

/ 2ixr 2 \ 2 

(/(r,q,t)/(r',q' l t)>= [yj^) ^, 9 o<W fe+fe+&6 > x 

e -ik ± (r ± -q ± ^i/g )-ik'_ L (r'_ L -q'_ L ^/9o) e -(fei+fef )r2/8-(g2 +^ 2 ) r 2 /2 ^ 

k ± ,k' ± 

where 6 + , 6 are the operators of the generated mode. The effect of the phase 
screen is represented in Eq. ( FlUj) by the parameter r\, determined via the 
correlation length, A c , of phase variation due to the phase screen as 

r 2 L2 fin 

In the absence of a phase screen, we may set formally A c = oo. Then it 
follows from Eq. (ITTj) that in this case r% = r . For any finite value of 
A c , ri < r and, as follows from Eq. ffTDl) . the characteristic values of the 
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transverse momenta of photons q±, q' ± are increased. This means that the 
beam becomes more divergent after passing a phase screen . 

With the known explicit terms Qj and ffTUj) . the calculation of (R^) 
reduces to many-fold integrations which can be performed straightforwardly. 
The result is given by 

(Rl) = 0.066n 2 rQcy^(Qon) 1/3 h. (12) 
The dimensionless quantity I\ is determined by the integral 

h = t dx(x - 1) V + (r 2 /4 + l' 2 ) g y iZ -r 1/6 , (13) 
Jo 

where 1' = Iq/(2tt). It can be calculated numerically. In the limiting cases of 
short and long propagation distances it is given by the following analytical 
expressions: 

h - ~ (l'o 2 + r 2 /4)- 1/6 , when + r 2 /4) » 1; 

o \q r i J 

27 

h « — , when qlr\z-\l% + r 2 /4) « 1. (14) 

Usually (?"q/4) >> Iq. Then, the above criteria mean that the diffraction 
broadening is smaller (upper case) and greater (lower case) than the initial 
beam radius. The upper case in Eqs. (IT4"|) results in (R 2 ) = 1.919C 2 2: 3 (2ro) -1 / 3 , 
that coincides exactly with the classic formula presented in the reference [TT] . 
(See Eq. (45) there.) As we see, there is no dependence of beam wandering 
on the phase screen when the propagation distance is very short. The result 
is evident for this limiting case in view of the fact that both the diffraction 
broadening and the broadening due to the atmosphere turbulence is much 
smaller than the initial radius of the beam. With the increase of propagation 
distance, z, or decrease of the initial coherence, the upper case in Eqs. (|T4l) 
may transform to the lower case which corresponds to dominating diffraction 
broadening of the beam. Then the dependence (R^) ~ r^ 3 will arise. As 
we see (R^) decreases with decreasing initial coherence. In this case the 
variance of the wander distance can be controlled by a suitable choice of the 
random phase screen. 

The situation is much more complex when the turbulence is strong. The 
averaging is no longer decoupled in the manner shown in Eq. ([6]). An essen- 
tial dependence of the distribution function on turbulence takes place here. 
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Therefore the approach based on employing Eq. (j5J) is not advantageous. 
The simplest way for further analysis is to proceed from the initial defini- 
tion of the wandering given by Eq. ([T|). The expression for the distribution 
function /(r, q, t) is given by [6J 

£1 jS 1 -ik ± jr-c q t+ — f ' dt't'Fj[r(t')l } ,+ , i 

/(r,q,t) = ^E e 1 90 J ° L J/fo S+k ± /2, go fo Q-k ± A g o/2| i =o, 

(15) 

where q = {q±,?o}j Q — q± — Jo ^'F±[r(t')], and r(t') is the trajectory of 
particle, which has the velocity c^u') an d is affected by the force F. The 
initial conditions are given by r(t' — t) — r and q(t' = t) — q. 
Substituting Eq. (TT5"1) into the general expression 



1 1 dr ± dr , ± r ± r^(/(r,i)/(r , ) t)), 



which determines the mean-square variation of the wander distance, and 
averaging over phase variations introduced by the random phase screen, we 
arrive at 



r' e~ ik±<r ~ Cq *^* k '- L ( r '~ c< i'*) 



x e -(fci+fcf)^/8 < / e -(Q 2 +Q' 2 W/2-(ic/g ) / t dt't'[k ± F(r(t'))+k' ± F(r'( i '))] \ ^g) 

where the last averaging is over random values of the refraction index. It is 
worth mentioning that momenta Q and Q' depend linearly on F. Therefore 
it is convenient to rewrite the quantity in the last angle bracket in a more 
convenient equivalent form as 



dpdpVto+p'^-fr +p/ )/(2r * } 



(27rr 2 ) 2 

; Jo d *'[(^ k ^ «'+p)F[r( i ')]+(^k' x t'+p')P[r' (*')]] \ 



x 

where p and p' are vectors perpendicular to the z-axis. 

Thus, the problem is reduced to the calculation of manyiple integral. 
There is a 13-fold integration in Eq. ([TBI) . After substituting Eq. (TTTj) into 
( fTBT) . the number of integrations increases to 17. Besides that, averaging over 
fluctuations of the refractive index introduces four additional integrations. 
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10 13 *C 2 (rrf 2/3 ) 



Figure 1: Dimensionless mean square wandering radius vs the turbulence 
strength for different values of the partial coherence determined by the 
ratio r\jr\. For all curves go = lC^m -1 ,^ = 10~ 3 m. The symbols show 
the results of numerical calculations. For visual convenience, the solid lines 
connect symbols using B-spline approximations. 



Finally, we have a 21-fold integral. We have performed most integrations 
analytically. The corresponding procedure is similar to that described in Ref. 
[6]. The rining, 5- fold integral, has been calculated numerically. Fig. 1 shows 
the results for fixed values of the aperture (r = 4) and propagation distance 
[z = 10km). We plot the dependence of the dimensionless quantity (R^) / R% 
on the turbulence strength. (It is the ratio (R^/Rf rather than merely (R^) 
that is informative about the practical significance of the wandering.) The 
beam radius Rb is given by the expression [6] 



c>2 _ r o 
tih — — 



4z 2 8z 3 T 



2 2 2 

Won 



(18) 



where T = 0.558C2/" 1 / 3 . 

As we see in Fig. 1, there is still considerable beam wander even for very 
strong turbulence, i.e. for C\ = 5 * 10~ 13 m -2 / 3 . (Usually, the value C\ = 
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10 _14 m _2//3 is considered to be a moderate turbulence level.) The four curves 
merge into a single curve when C 2 — > oo. This means that the wandering does 
not depend on the initial coherence in this case. So, the universal behavior of 
the wandering effect corresponds to the general concept that the atmosphere 
controls beam parameters for long-distance propagation or for the strong 
turbulence regime. At the same time, we see here the general tendency of 
the wander to decay with the increasing turbulence strength supports the 
reasonings of Fante [ID] . He considered that when the turbulence is strong, 
the beam no longer wanders significantly, but rather breaks up into multiple 
beams. 

In the opposite limiting case, C\ — > 0, the wander distance R w also tends 
to zero due to the obvious fact that the wandering is entirely caused by tur- 
bulence. From a formal point of view, there should be at least one maximum 
in the curve which connects the regions of weak and strong turbulence. The 
corresponding physical picture can be explained in terms of two competitive 
tendencies occurring when C 2 increases: (i) in the range of weak turbulence, 
where the beam radius is almost independent of the turbulence, the probabil- 
ity to meet sufficiently strong large-scale fluctuation of the refractive index, 
which deflects the beam as a whole, increases linearly with C 2 , (ii) in the 
range of strong turbulence, there is considerable beam widening due to pho- 
ton scattering on fluctuations of the refractive index (B% ~ C%) ; therefore 
the previous possibility has a low probability. This explains the presence of 
the maxima in Fig. 1. 

It is interesting to compare the results of the weak-turbulence theory given 
by Eq. ffl2l) with those of more general approach based on the distribution 
function (fT5]) . The results are shown in Fig. 2. As we see, both approaches 
give almost coinciding data for small values of C 2 . When C 2 increases, the 
results of the weak-turbulence theory are overstated. A similar picture was 
observed in Ref. [12] where the weak-turbulence theory was tested by means 
of computer simulations. 

Fig. 3 illustrates the dependence of beam wander on C\ for a shorter 
distance (5km) than in Fig. 1. The plots in Figs. 1 and 3 are very similar 
with the only difference that the maxima in Fig. 3 are displaced to the 
range of greater values of the turbulence strength C 2 . This difference is 
quite evident. Namely, initially the overall effect of the turbulence increases 
with the increase of both the value of C\ and the distance z. Therefore, 
the decrease of one of the factors can be compensated by the increase of the 
other one, thus providing almost the same effect of the turbulent atmosphere. 
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Figure 2: Dependence of beam wandering on the turbulence strength for 
small values of C\. Dashed lines show the dependence given by Eq. f[T2"j) , 
solid lines with symbols show the results of more general theory. 

Fig. 4 illustrates how the two approaches correspond to one another at small 
values of C\. Again we see a good agreement of both theories in this range 
of CI 

Fig. 5 illustrates the dependence of the ratio R 2 w jR\ on the turbulence 
strength for small values of the aperture radius (ro = 1cm). There is a 
significant decrease of the wandering effect in this case. This is because 
a small value of ro (and automatically ri) results in considerable diffraction 
broadening of the beam for such a long propagation path (5km). That is why 
the influence of turbulence on the beam parameters becomes competitive at 
greater values of C\ and, in correspondence with the latter, the maxima of 
both curves are displaced to the right as compared with Fig. 3. Also, the 
effect of partial coherence is more pronounced for smaller initial radius of the 
beam. This can be seen by comparing Figs. 3 and 5. 

The results presented in Figs. 1-5 require additional comments. Our anal- 
ysis proceeds from Eq. ([2]) where the evolution of the distribution function 
is based on the kinetic equation (TjJ. By definition this function is quadratic 
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Figure 3: The same as in Fig. 1 but for z = 5km. 

in the field amplitudes describing the intensity of the irradiance. Not all mo- 
menta and wave-vectors kj_ contribute to the observable intensity. One 
can easily see that initially in the absence of a phase screen, the characteristic 
values of both q± and k± are given by r^ 1 . In the presence of phase screen, 
the characteristic value of q± is of the order of rf 1 that determines the di- 
vergence of the beam and its broadening (diffraction broadening). It follows 
from geometric consideration that the diffraction broadening is of the order of 
z 2 5o~ 2 rJ~ 2 . This is almost the same value as given by Eq. (|18p . Therefore the 
characteristic value of k±_ decreases with distance as (r 2 /4 + z 2 rf 2 (fo~ 2 ) -1 / 2 . 
Also, the momentum of the moving particle varies with distance due to 
scattering on atmosphere inhomogeneities. The additional momentum ac- 
quired in this way, Aq^, can be estimated from its mean square value as in 
the case of a Brownian particle moving in q^-space and being affected by a 
random force Fj_ during the time t = z/c. Thus we have 

(Aql) ~ (Fl)t ~ Clz. 

As a result, the beam becomes more divergent and additional broadening 
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Figure 4: The same as in Fig. 2 but for z = 5km. 



due to the turbulence, Ai? 2 , can be estimated as 



ARl 



dz'Aq±(z') 



2\ 

C„z 



1 

Qo 

. This again agrees with Eq. (QJ), where Ai? 2 = 2.23/q 1/3 C^ 3 . 
When 

AT? 2 » | + ^rr 2 , (19) 

one can say that the beam size is determined almost entirely by the effects of 
the turbulence. In this case the characteristic values of k± are of the order of 
(Ai?fe) _1 and decrease with the increasing turbulence as C" 1 ^" 3 / 2 . Also the 
characteristic value of q± becomes of the order of Aq± which is much greater 
than its initial value r^ 1 (or rf 1 ). The last point can be seen directly from 
Eq. (TT9T) when we represent the turbulence broadening as 

ARl ~ Z - 2 Aq\. 
% 

The condition Ag_i_ >> r Q l -> r i 1 means a considerable randomization of the 
radiation field. The waves acquire properties of Gaussian statistics that is 
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Figure 5: The same as in Fig. 3 but for r = 1cm. 

very important when calculating beam wander variance. In contrast to cal- 
culations of the beam radius which is determined by correlations of only two 
waves, the beam wandering effect is determined by four-wave correlations (or 
by pair correlation function of the intensity (II))- The results presented in 
Figs. 1-5 were obtained explicitly assuming the dominating contribution to 
the average 

</(r)/(r')> -EE e-^ r+k ' r ')(6+ +k/2 6 q _ k/2 6j +k72 V- k 72) 

q,k q',k' 

to be from small regions of k and k' as explained above. (To simplify the 
notation, we omit the indices (j_) in all variables.) In this way strong cor- 
relations of pairs of waves &q +k / 2 ? &q-k/2 and frq/ +k // 2 , &q'-k'/2 were taken into 
account. At the same time it is evident that there is another region of wave 
vectors, i.e. 

|q + k/2 - q' + k'/2|, |q - k/2 - q' - k'/2| ~ AR^\ 

where pair correlations of waves may be also essential. The waves from differ- 
ent pairs, shown above, may correlate in this region. Conventionally, we will 
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refer to this type of correlations as cross-correlation. In the case of strong 
turbulence, the contribution of cross-correlations is not small, thus provid- 
ing saturation of fluctuations at high level. (See, for example, [6].) The two 
regions of wave vectors are well separated from one another and possible over- 
lapping in the course of summing over wave vectors is not important in the 
case of strong turbulence [T3]. When the turbulence effect becomes weaker, 
these regions approach each other, and in the limit of small turbulence they 
unite into a single region. In this case the beam wander is determined by the 
asymptotically exact solution (fl2l) . 

For strong turbulence, the contribution of cross correlations to the beam 
wandering, (i?^) cro ss) can be obtained as done in previous calculations. It is 
given by 



z 



2 



(R-w) cross „ 2 2 A r>? " (^0) 



w/cross 3r 2 q 2 AR 2 

Let us compare the value (R^) cross to (R^), shown in Figs. 1-5. First of 
all, consider those C 2 which correspond to the maxima in the curves plotted 
in Figs. 1,3,5. For the case r = r\ we see that AR 2 >> in all 

cases. This means that these maxima are in the range of strong turbulence, 
and Eq. fl20l) is applicable here. The values obtained from Eq. (120]) consist 
of only 7%, 5%, 0.4% of the corresponding data in Figs. 1,3,5, respectively. 
Moreover, if one moves towards greater values of C 2 , the contribution of 
cross-correlations will become of smaller because of increasing of AR 2 . A 
similar situation occurs when r\ becomes less than r$. 

On the other hand, our solutions with cross-correlations neglected almost 
coincide with those given by weak-turbulence theory when C 2 — > 0. (See Figs. 
2 and 3.) This assures us that Figs. 1,3,5 represent reasonable solutions for 
the specific set of parameters used there (and close to those) for any values 
of the turbulence strength, C 2 . 



3 Conclusion 

We have applied the method of photon distribution function [6\ to describe 
beam wander in turbulent atmosphere. In the limit of weak turbulence and 
in the absence of artificial random phase modulation, it becomes possible 
to obtain an analytical expression for the wandering radius, which coincides 
with the one known in the literature. Also, by bringing together analytical 
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and numerical calculations, we have succeeded in obtaining the wandering 
radius in the range of strong turbulence. The general conclusion of the actual 
studies is that the variation of the initial spatial and temporal coherence pro- 
vides significant positive (from the viewpoint of practical applications of laser 
beams) influence on the character of the intensity fluctuations. Namely, the 
relative value of the wandering radius can be considerably reduced. Moreover 
this reduction takes place just in the range of the most pronounced wander- 
ing effect. (See Figs. 1,3,5.) At the same time, the effect of partial coherence 
vanishes for very strong turbulence. This is in contrast to the behavior of 
the scintillation index, which in this case can be significantly suppressed by 
decreasing the initial coherence of the light. (See, for example, Refs. [I] and 
[6].) But this suppression is not very important because of small wandering 
effect in this case. 
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